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Nonlocal Gluon Condensate from the Dyson— Schwinger Equations 
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We establish a novel method of obtaining the gauge-invariant bilocal condensate of gluon fields 
from the solutions of the Dyson-Schwinger equations (DSE) which works in the infrared and ultra- 
violet limits. We present explicit results for 517(3) case in four dimensions, and compare to other 
related analytical and lattice results. 
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Motivation: The non-perturbative QCD vacuum is 
known for the presence of gluon (and quark) conden- 
sates These objects are in general studied in a lo- 
cal approximation. A bilocal object can be introduced 
by using operator product expansion (OPE) 2\. How- 
ever, OPE method can work only at a short scale. The 
infrared (IR) behaviour of gluon condensate is of impor- 
tance because it acts as a potential in which quarks fluc- 
tuate 0. Therefore, to understand, e.g. how QCD vac- 
uum properties are modified in the infrared by external 
electromagnetic fields one must first understand the 
dynamics of the gauge field condensate and than quark 
condensate in a scheme operational in the infrared. 

We obtain here the gauge-invariant nonlocal con- 
densate in terms of perturbative and non-perturbative 
propagators. Then we solve numerically the Dyson- 
Schwinger equations (DSE) for SU{3) gauge theory at 
approximately two loops, following the methods devel- 
oped by Bloch Q and other authors to obtain the non- 
perturbative propagators. We introduce perturbative 
propagators in a way which is consistent with the non- 
perturbative ones, and extract the gauge field conden- 
sate. Finally we analyze the condensate dependence on 
distance, and compare with prior results obtained by an- 
alytical and lattice methods. 

Form factors and nonlocal condensate: We start 
with the SU{Nc = 3) gauge field action 

I d'xG.^G'^-^, (1) 



where G^, = G^^T'^ , G^^ = d^Al 



d,Al+gr''-A''A\ 



. The gluon propagator Dp^^^ is consid- 
ered in Landau gauge Df^^{x) = -i(r (A°(x)A|^(0))) 
with a form factor F{p^) 
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and the ghost propagator is defined as D^°'^{x) 
-i{T (c'^(a;)c''(0))), with the ghost form factor G(p2), 

gab 



(2) 



-Gip')- 



(3) 



Following Ref. [6[, the bilocal condensate is defined 
as the vacuum expectation value (VEV) of the following 
gauge-invariant operator 



a tr 



[G^MW{y,z)G^^{z)W{z,v)] :), 



(4) 

where Wilson-Polyakov lines are W{y, z) = 
Vexp {ig Jy A^dx'^^ , V being the path-ordering symbol. 
Note that the normal ordering in Eq. ^ is defined with 
respect to the perturbative (PT) vacuum, whereas the 
operators are averaged over a non-perturbative (NP) 
vacuum. 

We are going to work at below the two- loop limit, thus 
Eq. ([4]) is expanded perturbatively up to a^-order: 
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indices C^),'-'^' denoting a order. The two-gluon piece is 
^^'^ = (■■ [GL'"^'^(2;)G;;jf (0)] :), (6) 

where we have split the gauge covariant quantity G^^ 
into a linear (Abelian) and a bilinear non-gauge covariant 
objects Gj^Jj^ = G^„ + ig[A^,A^]. Three-gluon term is 
absent, and a four-gluon term without Wilson lines is 



(1) _ 



a tr 



: - {g[A,{x),A^{xMA,{0),A,m} :). (7) 



Wilson lines do not non-trivially contribute to the above 
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and but contribute to J^j^-* and .7^3^'' ■ 



'F^,{x)Af3 {u{t))A^ (via)) x 



X F^"'{0)u'^v'^9{T-cr)dTda 
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a tr 



g^F^,^{x)Asi (u(t)) x 



X F^'''{^)A^{v{G))iPvPdTda 



(8) 



2 



We use here straight Unes v'' {a) — ax'' , (r) = tx^ , 
running over re (0, 1), (t G (0, 1). 

To express !F^^^ , T^l 3 in terms of propagators, we need 
to handle the PT-normal products in the NP-vacuum. 
Consider an operator 0{x). A T-product of opera- 
tors can be evaluated in both perturbative and non- 



perturbative vacuum, with Gq ' defined as follows: 



T{0{x)0{y)) = : 0{x)0{y) +iGPJ{x,y) 
= ■.0{x)0{y) -.^P +iG^P{x,y), 



(9) 



indices NP and PT over normal ordering symbols and 
Green's functions meaning the definitions of the latter in 
the nonperturbative or perturbative vacuum correspond- 
ingly. Taking a VEV over the nonperturbative vacuum, 
and subtracting the NP definition of T [0{x)0{y)] from 
the PT one, one sees that 

(: 0{x)0{y) ■.p^)mp = *(Gg^(a;, y) ~ Gg^(x, y)). (10) 

Thus emerges the object of interest 

^ o.[j?) {F^^ip') - F^^ip')) . (11) 

In our case the non-perturbative form factor 
pNP^p2^ _ pij)"^^^ Eq, is what we obtain from solv- 
ing the DSE. The corresponding perturbative form factor 
pPT (jp.-^ is discussed below. For the sake of generality 
we present, using Eq. the components of condensates 
in an arbitrary dimension d: 



=2((i-l)4 
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'^Apip')e^P^, (12) 



in agreement with Ref. 0], where a local limit of this 
equation is given, here and below dc = — 1. For the 
subleading terms we obtain: 

.ff ) = M,N, [{d^ -2d- l)/2 - , (13) 

where the propagator moments are defined as 

hix) = /^W),.p. 



For the terms containing Wilson lines, one has 

the moments Wi , W2 being 
d'^p 



Wi 



-A^^(p^)e'f^ 



(14) 



(15) 



(16) 



W2 = 12(1 -d) 



d-^p Apip^)f 2 (P2:)2\sin2(f ^ 



To be able to complete the evaluation of the nonlocal 
condensate we should supply the non-perturbative and 
perturbative form factors into these equations. 

Dyson— Schwinger equations: To obtain the needed 
non-perturbative input form factor F^P{p'^) one option 
is to solve DSE. We work in pure gluodynamics, following 
closely Bloch's method jSj of IR analysis and his trunca- 
tion labeled (A.l), Table 3 ibid., thus we do not give 
here the full details. The procedure of writing down, 
truncating, IR-analyzing and solving DSE is developed 
and described in detail in many works, see for exam- 
ple [i, B E, [m, [i3- The DSE equations in a covariant 
gauge are organized in the following way: 



4^-^ = -(n(^,2)-n(.|.) 



(17) 



where Ug and Uc are renormalization points for gluon 
and ghost propagator respectively. After truncation, the 
ghost self-energy and the gluon vacuum polarization are, 



E(p2) = N,gl^T,{p\q\r^)G{q^)F{r^)^ 



(18) 



where = yj Aira^jj), ji is renormalization point, and 
separate contributions to vacuum polarization are 

n''(p') = /(£f.A/o(p^g^r2)G(q2)G(r2), 

^''ip') = I^Qo{p',q^r')G{q^)G{r% ^^^^ 

n^^ip') = /(£f.g4(p^9^r2)x 

xG{qyF{qy/^G{r^)F{rY^\ 

with kernels Mo,Qo,Q4i?o defined in i, eqs. (47, 48, 
49, E.l). Here x = p^ ,y = q^ , z = r^, r'^ = (p — q)^. 

We decompose the gluon and ghost form factors in a set 
of Chebyshev polynomials Ti (z) on an interval p^ G (e, t) , 
and use a power-like Ansatz for the IR domain p^ E (0, e): 



F{z)^ 



'i^PiT,ifiTi{z)),z > e 
Az^'^, z < e 



(20) 



and a similar decomposition is assumed for G, where 

2 

z — logT-, uo is scale unit, continuity is imposed at 

Mo 

p2 — e, renormalization points are chosen for gluon and 
ghost respectively, af — a,(Tg — Q, normalization condi- 
tions F{(j) — 1, 1/G(0) = 0; and from now on, d = A. 
For the two sets of coefficients we get algebraic equa- 
tions, which are solved by Newton's method. Reliability 
of the method is supported by convergence of the sets of 
coefficients. 

Choice of scale: Upon solving the equations, we con- 
struct running coupling out of the form factors from DSE 



(pxf 



a{p')^a{ti)F{p^)G\p'). 



(21) 
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This follows since we did not include the ghost-ghost- 
gluon vertex factor into the system of DSE, it is consid- 
ered to be unity. Here fi is the rcnormalization scale for 
coupling, which is identical with the ultraviolet (UV) 
cutoff in our case: ^ — cr, and F{^^) — G{p?) = 1 

To relate the gauge coupling constant a with physics, 
one must introduce a scale. Physical scale is not fully 
defined in this problem, for we cannot actually compare 
pure gluodynamics to real-world QCD with fermions. 
Even so, we note that any object of dimension k must 
be measured in units of /Iq, where scale /^o is one and 
the same for the whole problem. By comparing our 
(non-perturbative) a with the standard running coupling 
value [l3| we obtain a rough estimate for our scale 



a{-s/s = lOGeV) = 0.18, ^lo ^ 0.55GeV. 



(22) 



For comparison, an estimate performed at ^/s — rrir = 
1.77GeV, a = 0.35 yields fiQ = 0.8. 

From DSE to condensates: Our solutions of DSE 
confirm the results by Bloch in every feature; thus we 
proceed to applying them for obtaining the conden- 
sate. The form factors which we have obtained and 
the running coupling are true non-perturbative quantities 
F^P{z),G^^{z),a{z). Now, in order to employ Eq. d!]), 
we need the difference between the perturbative and non- 
perturbative objects. 

Bloch in Q observes that he reproduced the leading- 
order PT results for coupling in the UV, but not the 
NLO. We confirm this by comparison of f3{a) and anoma- 
lous dimension 7(0;) for the propagator, see below, with 
their perturbative behaviour. We believe this is a conse- 
quence of the partial rather than full resummation of the 
perturbation series, since not all of the NLO diagrams 
emergent from the systematic expansion appear in DSE 
as solved here, due to the truncation of the hierarchy of 
DSE. In this sense we might call the present approxima- 
tion resummed 1.5 loops. Resummed, since the partial 
resummation extends to all orders. For possible improve- 
ments of our method, it would be advantageous to use 
DSE with vertex function included [l^ . 

Given that the DSE solution is of infinite order in a, 
but incomplete as of 2nd order, one must re-evaluate per- 
turbation theory form factor, rendering it consistent with 
the partial resummation inherent in the truncated DSE. 
This must be always done since DSE is always in practical 
work a truncated system. Gluon form factor, normalized 
at fj,, is expected to be organized in the form: 



F(p2)^exp / 



b{a) 



da 



(23) 



where h{a{z)) — a{z)^{a{z)) / (3{a{z)), z, as before, is 
renormgroup time. Functions /3(a) and 7(a) are nonper- 
turbative and given by: 



P{a{z)) 



da{z) 
dz 



7(a(z)) 



dlog F{z) 
dz 



(24) 
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FIG. 1: The nonlocal gluon condensate in 3 + 1-dimensional 
517(3) from DSE approach, scale fixed at 10 GeV. 



In perturbation theory IJ] b^^{a) — ^ J2i=o ^i^^j where 
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ever instead of these hf^ we must use results of a poly- 
nomial fit in the UV range of our ratio b{a) 



a7(a) 



b{a) ^Yka\z) 



PT 



i=0 



(a). 



(25) 



To obtain a reasonable precision, n ~ 2 suffices. Using 
/3, 7 from DSE and having performed the fit Eq. (|25|) . we 
apply the coefhcients bi to build an "effective" pertur- 
bative expression for gluon form factor, which is to be 
afterwards subtracted from the NP expression: 



F^^(p2) = exp 



(26) 



We have now got both, a non-perturbative propaga- 
tor from solving DSE, and the perturbative propagator, 
Eq. ([26| . Thus we can obtain the nonlocal condensate 
shown in figure [1] This is our main computational re- 
sult. Our solution qualitatively satisfies the expectations 
which are due to the known UV and IR forms of the 
condensate, prescribed by OPE and lattice respectively. 
OPE predicts at small distances [§| 



(27) 



Comparison with this short-distance behaviour is shown 
in table |T1 At large distances one can use [15| 



Tlx) = Ce- 



(28) 



(to be used with caution, as some other predictions e.g. 
from instanton models [l6|, exist). It is instructive to in- 
spect numerical convergence to the values of G§, A^, C, a 
shown in table [III £^s function of the number n of poly- 
nomials involved. Results are stated for scale fixed at 10 
GeV. The remaining numerical error in the determination 
of the magnitude of condensate and correlation lengths 
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TABLE I: Comparison our results with prior work. "SR" 
denotes sum rules, r means "r decay" . Results of our present 
Letter given in bold in the two last lines, denoted by f and 
X, for scale fixing points 1.8 GeV and 10 GeV respectively. 
Error bars arc given in parentheses. 



Method 


Yr. 


Ref. 


Gg[GeV-] 


A^[GeV"] 


C[GeV*] 


a [GeV] 


SR 


78 




.012 








OPE 


'92 


m 




.21 






Lattice 


'97 


\17] 


.015 




.008 


.6 


Lattice 


'02 


\18] 






.08 


.8 


SR 


'02 


\19] 


.009(.007) 








T 


'02 


[19] 


.006(.012) 








DSE 


'08 


t 


.062(.022) 


2.3(.6) 


.09(.03) 


1.9(.4) 


DSE 


'08 


t 


.014(.005) 


1.1(.3) 


.020(.007) 


1.3(.3) 



Jacques Bloch for introducing us into the details of 
his DSE method. We thank A. Bakulev, A. Dorokhov 
and S. Mikhailov for discussions. Numerical calculations 
were performed in part on the Computational Cluster 
of the Theoretical Division of INR RAS, Moscow. A.Z. 
would like to thank D.V. Shirkov and A.S. Gorsky for 
their kind advice. We thank P.D. Dr. Peter Thirolf and 
Prof. Dr. D. Habs, Director of the Cluster of Excel- 
lence in Laser Physics - Munich-Center for Advanced 
Photonics (MAP) for their hospitality in Garching while 
much of this research was carried out. This research 
was supported in part by RFBR Grant 07-01-00526, by 
the DFG-LMUexcellent grant, and by a grant from: the 
U.S. Department of Energy DE-FG02-04ER4131. 



TABLE XL Convergence of condensates and correlation 
lengths with respect to the number of polynomials n. 



n 


G;^[GcV-] 


A^[GcV"] 




a[GeV] 


22 


0.00834 


1.04052 


0.02671 


2.26331 


24 


0.01068 


0.64805 


0.01854 


1.22284 


26 


0.00967 


1.30711 


0.01573 


1.41736 


28 


0.01270 


1.18162 


0.01892 


1.33049 


30 


0.01352 


1.12676 


0.01977 


1.28204 



can be eliminated by increasing the size of the basis, once 
there is merit in such an effort. Bloch used n = 50 
in solving DSE. We estimate our current numerical error 
based on convergence shown in table [IT] and show this in 
the bottom entries in table HI 

Conclusion: We described here what we believe is a 
novel idea allowing to obtain the non pcrturbative vac- 
uum condensate, which has the capability of yielding 
both IR and UV vacuum properties. This approach is 
less expensive than lattice simulations, and more univer- 
sal than OPE. 

Our main result is the demonstration that the solution 
of DSE allows determination of the full nonlocal gluon 
condensate (: ^§-G ^^{x)W{x,G)G ^^{0)W{0,x) :) as a 
function of coordinate or momentum. From these one 
can get its fundamental properties: the value at zero sep- 
aration and correlation range. Our results are obtained 
in Landau gauge, yet the procedure and thus the final 
answer are gauge-invariant. Our here presented method 
and results suggest that the effort required, especially in- 
cluding Fermi (quark) fields, is greatly reduced compared 
to lattice method. 

Moreover, in principle there seems to be no obstacle to 
the study of the far infrared limit, which relies of further 
development of DSE in this domain. Our pure gauge the- 
ory results cannot be as yet used for further phenomeno- 
logical developments due to scale uncertainty; however, 
the purpose of this paper was to demonstrate in principle 
how an a priori calculation of condensate is possible in 
the framework of DSE. 
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